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Feasibility study of using the overlap-Dirac operator for hadron 
spectroscopy* 

UKQCD Collaboration, Craig McNeile, Alan Irving, and Chris Michael, a 
a Department of Mathematical Sciences, University of Liverpool, L69 3BX, UK 

We investigate a number of algorithms that calculate the quark propagators for the overlap-Dirac fermion 
operator. The QCD simulations were performed at f3 = 5.9 with a lattice volume of 16 3 32. 



o 
O 

(N 

(N 
> 
OS 

in 
o 

Os 
O 
Os 
OS 
-t— > 

Oh! 

<D 

43 



X 



1. INTRODUCTION 

Two critical systematic errors in the calcula- 
tion of the fs decay constant from lattice QCD 
are the chiral extrapolations and the unquenching 
errors jjj. The only way to reduce these errors 
is to simulate QCD with lighter quark masses. 
Unfortunately, because of exceptional configura- 
tions, it is difficult to further reduce the masses 
of the light quarks in quenched simulations with 
the clover operator. Progress in reducing the sea 
quark masses in dynamical fermion simulations 
with Wilson like quarks is slow |Q. 

It seems plausible that the difficulty of simulat- 
ing with light quark masses with the clover oper- 
ator is due to explicit chiral symmetry breaking 
in the action. Neuberger has derived a fermion 
operator, called the overlap-Dirac operator, that 
has a lattice chiral symmetry . 

Our goal is to simulate the overlap-Dirac oper- 
ator in a mass region: (Mps/My = 0.3 — 0.5), 
which is inaccessible to clover quarks (but not 
staggered quarks J| ) . Most of the techniques de- 
veloped in the quenched theory can be used for 
full QCD simulations B. 



2. THE OVERLAP-DIRAC OPERATOR 

The massive overlap-Dirac operator HQ is 



„ m 1 / , n Him) , . . 

D N = -(1 + M + (1 - M)75 \l ) (1) 



where H(m) is the hermitian Wilson fermion op- 
erator with negative mass, defined by 

H(m) — J5(D W — m) (2) 

where D w is the standard Wilson fermion oper- 
ator. The parameter fj, is related to the physical 
quark mass and lies in the range to 1. The m 
parameter is a regulating mass, in the range be- 
tween a critical value and 2. The physics should 
be independent of the mass m, but the value of m 
effects the locality of the operator and the number 
of iterations required in some of the algorithms 
used to compute the overlap-Dirac operator. 

3. NUMERICAL TECHNIQUES 

Quark propagators are calculated using a 
sparse matrix inversion algorithm. The inner step 
of the inverter is the application of the fermion 
matrix to a vector. For computations that use 
the overlap-Dirac, the step function 

e(H)b = -£=& (3) 

must be computed using some sparse matrix al- 
gorithm. The nested nature of the algorithm 
required to calculate the quark propagators for 
the overlap-Dirac operator makes the simulations 
considerably more expensive than those that use 
traditional fermion operators. 

Practical calculations of the overlap operator 
are necessarily approximate. To judge the accu- 
racy of our approximate calculation we used the 
(GW) Ginsparg- Wilson error: 
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which just checks that the matrix obeys the 
Ginsparg- Wilson relation pj. 

Our numerical simulations were done using (3 = 
5.9 quenched gauge configurations, with a volume 
of 16 3 32. The quark propagators were generated 
from point sources. For all the algorithms we in- 
vestigated, we used m equal to 1.5. 

4. LANCZOS BASED METHOD 

Borici || has developed a method to calcu- 
late the action of the overlap-Dirac operator on 
a vector, using the Lanczos algorithm. In exact 
arithmetic, the Lanczos algorithm generates an 
orthonormal set of vectors that tridiagonalises the 
matrix. 

HQ n = Q n T n (5) 

where T n is a tridiagonal matrix. The columns of 
Q n contain the Lanczos vectors. 

The "trick", to evaluate the step function 
(Eq. ||), is to set the target vector b, as the first 
vector in the Lanczos sequence. An arbitrary 
function / of the matrix H acting on a vector 
is constructed using 

(f(H)b)i = £(Q„/(T„)Qt). , & . (6) 
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where the orthogonality of the Lanczos vectors 
has been used. The f(T n ) matrix is computed 
using standard dense linear algebra routines. For 
the step function the eigenvalues of T n are re- 
placed by their moduli. Eq. [?] is linear in the 
Lanczos vectors and thus can be computed in two 
passes. 

The major problem with the Lanczos procedure 
is the loss of the orthogonality of the sequence of 
vectors due to rounding errors. It is not clear 
how this lack of orthogonality effects the final re- 
sults. Some theoretical analysis has been done on 
this method || . It is claimed that the lack of or- 
thogonality is not important for some classes of 
functions. 

In Fig. [I], we plot the eigenvalues of the 
Ginsparg- Wilson operator, as a function of the 
number of Lanczos steps, for a 2 4 hot 51/(3) con- 
figuration. As the number of Lanczos steps in- 
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Figure 1. Eigenvalue spectrum of the overlap- 
Dirac operator, for 100 (crosses), 50 (diamonds), 
and 10 (bursts) Lanczos iterations. 

creases, the eigenvalue spectrum moves closer to 
a circle (the correct result). Even after 50 iter- 
ations of the Lanczos algorithm, there are still 
small deviations from the circle. Unfortunately, 
it is much harder to look at the eigenvalue spec- 
trum using a production gauge configuration, so 
we computed the GW error instead. The GW 
error was: 5 10" 3 (50 iterations), 6 10~ 4 (100) 
iterations, and 3 10~ 4 (300 iterations) on a single 
gauge configuration. 

Fig. ||, is an effective mass plot of the pion, 
for two choices of mass and number of Lanczos 
sweeps. The "plateau" in the pion effective mass 
plot for approximate operator that used 50 Lanc- 
zos iterations is higher than the lowest pion mass 
that can be reached with non-perturbatively im- 
proved clover. It is not clear what causes the 
"shoulder" in Fig. |[ We would like to compute 
the eigenvalues of the overlap-Dirac operator on 
the bigger gauge configuration, to check how ac- 
curately we are computing the overlap-Dirac op- 
erator. 
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Figure 2. Pion effective mass plot 

5. RATIONAL APPROXIMATION 

The step function can be be approximated by 
a rational approximation. 



e{H) ~ H (co 
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The rational approximation typically approxi- 
mates the step function, between two values. The 
eigenvalues of the matrix H should lie in the re- 
gion where the approximation is good. The coeffi- 
cients Ck and dk can be obtained from the Rcmcz 
algorithm ]l0| . The number of iterations required 
in the inverter is controlled by the smallest dk co- 
efficient, which acts like a mass. We have not 
yet implemented the technique of projecting out 
some of the low lying eigenmodes . 

On one configuration we obtained GW errors 
of: 1 1(T 4 , and 4 1CT 5 , for the N = 6, and N = 8, 
optimal rational approximations p0| . The mul- 
tiplicative scaling factor was tuned to obtain the 
best results. Unfortunately, the above results re- 
quired up to 600 iterations for the smallest dk, 
which was too large to use as the inner step of a 
quark propagator inverter. 

One feature of the optimal rational approxima- 
tion JTofl , is that the lowest dk factor is smaller 
than the square of the validity of the approxima- 
tion, which means that the condition number of 
the inversion is that of the matrix H 2 . We ex- 



perimented with a hybrid quadrature and series 
approximation to Robertson's integral represen- 
tation of the step function. 



e(H) = 
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In Eq. ^, the first integral was approximated using 
an open quadrature rule and the second integral 
was approximated by a Chebyshev series. The 
step size in the quadrature formulae reduces the 
condition number of the required inversion. In 
our preliminary tests of the algorithm, the hybrid 
method produced a substantial reduction in the 
number of iterations required over the optimal 
rational approximation. However the computed 
solution was less accurate than that produced by 
the optimal rational approximation, because the 
Ginsparg- Wilson error was 6 10 -4 . 

Clearly more work is required on the algo- 
rithms that calculate the step function, before the 
overlap-Dirac operator can be used in the quark 
mass region we are interested in. 

This work is supported by PPARC. The com- 
putations were carried out on the T3D and T3E 
at EPCC in Edinburgh. 
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